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Abstract 

The paper considers some typical problems for kinetic models evolving through pair-collisions at 
temperatures not far from absolute zero, which illustrate specifically quantum behaviours. Based on 
these - mostly recent - examples, a number of differences between quantum and classical Boltzmann 
theory is then discussed in more general terms. 



Introduction. 

Let us start by recalling some basic facts about the classical Boltzmann case. Consider a parcel 
of gas evolving as 

dx dv 

-r = v, -r = F- 
dt dt 

That implies the evolution of the gas density /; 

n tt u\ n\ *\ d f , d f dx , d f dv d f , d f , j? d f 
D t f( x (t),v(t),t) = _ + __ + __ = _ + «_+ ir_, 

which gives the classical BE when driven by a collision term Q(f), 
(D t f(x(t),v(t),t) + v ^L + F^- = Q(f). 



For two particles having pre-collisional velocities u, «*, denote the velocities after collision by v',v^, 
and write /* = /(v*), /' = f(v'), f* = f{ v i) • I n this notation the classical Boltzmann collision 
operator becomes Q(f) = J B(f'f^ — f f*)dv*duj. The kernel B is typically of the type b[oS)\v — v*\ s 
with —3 < s < 1. The collision term is linearly proportional to the density of each of the two 
participating molecules. Mass and kinetic energy are each conserved, and the monotone in time 
entropy J dvdxf log f(t) prevents strong concentrations to form. In the classical case all collisions 
respecting the conservation laws are permitted, 
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The quantum case is different. Close to absolute zero, permitted energy levels are often discrete, 
implying fewer collisions. There may exist a condensate, where an excitation may interact only 
with those modes of its zero-point motion that will not give away energy. In the quantum regime, 
also the thermal de Broglie wave length may become much larger than the typical inter-particle 
distance, a situation totally different from classical kinetic theory. The scattering is wave-like and 
two-body quantum statistics gives a marked imprint on the individual collision process, e.g. the 
Pauli principle for fermions. 

The rich phenomenology of low temperature gases corresponds to a similarly rich mathematical 
structure for its models, including the kinetic ones. The quantum kinetic collisions and equations 
are much more varied and fine-tuned than in the classical case. This paper will consider typical 
quantum kinetic models evolving through pair collisions, in order to illustrate the differences be- 
tween the classical and quantum cases; 

1) The Nordheim Uehling Uhlenbeck model, 

2) Spin in the fermionic case, 

3) Low temperature bosons and a condensate, 

4) A kinetic anyon equation. 

These models are obtained by physics arguments from wave mechanics in the Heisenberg setting 
and from quantum field theory. In some cases the models have been validated in a formal mathe- 
matical sense (see [BCEP], [SI]). Based on the examples 1-4, the final section will elaborate on some 

5) Differences between quantum and classical Boltzmann theory. 



1 The Nordheim Uehling Uhlenbeck model. 

Around 1930 the first quantum kinetic equation for boson and fermion gases was proposed indepen- 
dently by Nordheim ([N]) and by Uehling-Uhlenbeck ([UU]). The equation is of particular interest in 
a neighbourhood of and below the transition temperature T c , where a condensate first appears. The 
NUU equation is semiclassical; the transport left hand side of this evolution equation is classical, 
only the collisional right hand side is directly influenced by the quantum effects. The term quantum 
Boltzmann equation usually refers to a semiclassical equation, even though true quantum kinetic 
models are sometimes considered. 

The collision operator for the NUU equation is 

QnuuUXp) = I B5(p + p* - (j/+£))8(E(p) + E{p*) - (E(p') + £7(pJ,))) 

JK 3 xR 3 xR 3 (11) 

(f'fl(l + e/)(l + e/„) - //*(1 + e/')(l + efl))dp*dp' 'dp'„ 

with e = 1 for bosons, e = — 1 for fermions, and e = gives the classical case. Here the quartic 
terms vanish. A discussion of the validation can be found in [BCEP]. 
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For e = ±1 the entropy is 

j (/log/-i(l + e/)log(l + e/))dp. 

Also this entropy is monotone but, contrary to the classical case, gives no control of concentrations. 
In equilibrium Q is zero, and multiplication with log and integration gives 

= / Q(f) log —j— f dv = [ B5(p + p* - (p'+p'*))6(E(p) + E(p*) - (E(p') + E&))) 

(1 + £ /)(l + e/,)(l + e/,(l + t /:)( 1+ / ' e/ , 1+ V r ^ r ^) 

The two last factors in the integral have the same sign, which implies 

/ /* / /* 

1 + ef 1 + eft = 1 + e/l + e// 

As in the classical case, we conclude that = M, a Maxwellian, hence / = , which for 
e = ±1 is the Planckian equilibrium function. 

For fermions (e = —1) the concentrations are bounded by one, because of the factor (1 — /) which 
preserves positivity together with f. This is stronger than the entropy controlled concentrations from 
the classical case. There are general existence results for the fermion case in an L 1 n L°°-setting 
due to J. Dolbeault [D] and P.L. Lions [PLL]. The boson case is more intricate. So far there are 
only partial existence results - for the space- homogeneous, isotropic setting ([Lul], [Lu2], [EMV]), 
and for some space-dependent settings close to equilibrium ([R]). 

Mathematical results have also been obtained concerning existence, uniqueness and asymptotic 
behaviour ([EM]) for the related Boltzmann-Compton equation 

k 2 ^ = J~(f'(l + f)B(k', k; 9) - /(l + f')B(k, k'; 0))dk> (1.2) 

with / photon density, k energy, 9 temperature, and with the detailed balance law e k ^ e B(k', k; 6) = 
e k / 9 B(k, k'; 9). This equation models a space-homogeneous isotropic photon gas in interaction with 
a low temperature gas of electrons in equilibrium, having a Maxwellian distribution of velocities. 



2 a) Spin in the fermionic case. 

The experimental study of spin polarized neutral gases at low temperatures and their kinetic mod- 
elling is well established in physics, an early mathematical physics text in the area being [S]. The 
first experiments concerned very dilute solutions of 3 He in superfluid 4 He with - in comparison 
with classical Boltzmann gases - interesting new properties such as spin waves (see [NTLCL]). The 
experimentalists later turned to laser-trapped low-temperature gases (see [JR]). The mathematical 
study of these models, however, is less advanced. To discuss this, we first recall some properties of 
the Pauli spin matrices 

CTl= (i o)' a2= (-i o)' CT3= (o 
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With [<7j,<Tj] denoting the commutator a^aj — (Jj(Ji, the Pauli matrices satisfy 

[oi,<72] = 2w 3 , [cr 2 ,cr 3 ] = 2ic7i, [cr 3 , cri] = 2icr 2 , and [crj,(Jj]=0 for i = 1,2,3, 

With cr = (<7i, o"2, (73) the Pauli spin matrix vector, this is equivalent to a x <r = 2icr. Let A4 2 (C) 
denote the space of 2 x 2 complex matrices, and H 2 (C) the subspace of hermitean matrices. % 2 (C) 
is linearly isomorphic to IR 4 , if we use the decomposition p = A C I + A s • a and identify p G "H 2 (C) 
with (j4 c , A,) G M 4 . 

A dilute spin polarized gas with spin ^, will be modelled by a distribution function matrix p G "H 2 (C) 
which is the Wigner transform of the one-atom density operator for the system being modelled. The 
domain of p(t,x,p) is taken as positive time t, p G M 3 , and for simplicity here with position-space 
x periodic in 3d with period one. We shall focus on the following simplified model for the kinetic 
evolution of p, 

d 

D P-= g- t P + P-VxP = Q(p) (2.1) 
with (in the Born approximation), the collision integral given by 

Q(p) = J dp 2 dp' l dp' 2 B5( Pl +p 2 - p[ - p' 2 )8(pl +p 2 2 - p'l - p'i) 

{[pi',/5i]+Tr(p 2 P2') - [pv , pi]+Tr(p 2 P2')} ~ g{[[Pi> Pi']+> \fc, P2'}+}+ ~ [[pi, Pv]+, [p2, p2'}+]+} 

Here p = I — p, and [.,.]+ denotes an anti-commutator. The kernel B is of hard force type, 
< B(\pi — p 2 |,#) = \pi — P2\^b{6) with < (3 < 1, C ^P^ G L°°. The number density of particles 
of any spin component, is given by / := Tr(p(t, x,p)), and the magnetization of particles is given 
by the vector a(t, x,p) := Tr(ap(t, x,p)), which implies p = \{fl + a ■ a). The resulting equations 
for / and a are 

Df = Q n (f,a), (2.2) 
Da = Q m {f,a), (2.3) 

where 

Qn(f, °) = \J dp2dp' 1 dp / 2 B5(p 1 +p 2 - P [ - P' 2 )t(pl +P 2 2~ P'l - P'i) 
(2 ([-ft' ~ 2^ 1 ^ 1 ' + dl ' ^ 1 ')][/2' - \U2f2' + 02 ■ G2>)\ ~ [fl ~ ^(flfl' + Ol ■ V )\\f2 ~ ^{hh' + °2 ■ OV)]) 
-^([^1' - ^(/l'^l + Ml')] • [°V - ^{h'02 + M2')] ~ [o-l - ^(/l'^l + Ml')]^ - ^{h'P2 + M2')])) » 

Qm(/, °) = \ j dp2dp 1 dp 2 B5{pi +P2- p[ - p' 2 )S(pl +pl~ p'l - p'l) 

(J([/2 - ^(/2/2' + 02 • ^ 2 ')][ai - + Ml')] - [/2' - ^Cfe/tf + ^ 2 ■ U 2 ')][a V - ^(/i'^i + Ml')]) 

~^([/l' - ^(/l/l' + ^1 • Ol')][oV - ^(/2'^ 2 + M2')] - [/l - ^(/l/l' +Ol • ^l')][^2 - ^(/a'^ + My)]))' 

The collision term Q n does not change the number density (J Q n dp = 0), the linear momentum 
density (J* pQ n dp = 0), and the energy density (J p 2 Q n dp = 0), and the collision term Q m does not 
change the magnetization density (f Q m dp = 0). 
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2 b) Existence results. 

Consider the initial value problem for the equations (2.2-3). The initial value (/o,^o) £ L°° is 
assumed for each ball p 2 < j 2 to satisfy < r]j < fo < 2 — rjj, min((2 — /o) 2 , /q) > oo • <7o + rj 2 for 
some rjj > 0. We first treat the case with a truncation Bj in the domain of integration for Q, where 
Bj is the restriction of B to the set p\ + p 2 , < j 2 . 
Set 

F(t,x,p) = f(t,x,p) for 0</<2, =0 for / < 0, = 2 for / > 2, 
£(i, x, p) = a(t, x, p) when min(F 2 , (2 — F) 2 ) > a ■ a, 

. min(F,2-FW, 

else ij{t,x,p) = (t,x,p). 

V<r ■ cr 

The system £)/ = Q n (F, £), Dct = Q m (F, S) with the initial value (/o, ^o) under the truncation I?,, 
can be solved by a contraction argument. 

Lemma 2.1 [A2] T/ie initial value problem with initial values (fo,o~o) for the truncated problem 
Df = Q n (F, £), Da = Q m (F, S) with truncation Bj, has a unique local solution in L°°. 

It remains to prove that F = f and X = a. This holds by continuity on a (short, n- and j-dependent) 
time interval < t < tj, using the boundedness of T n and T m . We give the proof for general t, 
when there is spin only in the 03 direction, and refer to the paper [A2] for the general case. With 
a = (0,0, s), we can replace a with s in Q n and Q m , and the equation for a by the corresponding 
one for s. Set c\ = 18max Pl f dp2dp' l dp' 2 B n 5{pi + P2 — v'\ ~ ^2)^(^1 +P2 — ~P'i)- ^ holds 

- - s 2 ) < j dp 2 dp' 1 dp' 2 B6( Pl +p 2 - pi - p' 2 )% 2 + p 2 - K 2 - p 2 2 ) 

^( " (/l " *l)([(/l + -l)(l " ^(/l' + -l'))][(/2 - «2)(1 " ^(/2' " **))]) 

+ *i)([(/i - -0(1 - ^(/r - «i0)][(/2 + «2)(1 - ^/ 2 < + **))])) 

+ ^ ( " (/l " «l)([(/l + -l)(l " + *l0)][(/2 + -2)(1 " \{f V + S2>))\) 

-(/l + *l)([(/l " -1); (1 " \(h ~ sy))][(f2 ~ S 2 )(l ~ ~ S 2 ,))])) 

< (fl ~ Sl)(Qn + Qm)(f, S) + (/l + ai)(g„ - Q m )(f, s) = D(f 2 - S 2 ). (2.4) 

This implies 

< (f 2 - s 2 )(x,p)e-^ < (f 2 - s 2 )#(t, x,p), (2.5) 

where (/, s) is the contraction solution with initial value (/o,<7o) • By continuity / ± s > for 
< t < tj. Analogously, starting from the equation for (1 — \{fi ± o"i)) instead of /i ± a\, we get 
by uniqueness the same solution, together with the estimate 

< ((2 - /o) 2 - s 2 Q )(x,p)e-^ < ((2 - f) 2 - s 2 )#(t, x,p), (2.6) 

which by continuity holds for < t < tj together with 2 — / ± s > for < t < tj. And so 
uniformly in (x,p), \s\ < min(/, 2 — /) < 1 for < t < tj. By iteration, existence and uniqueness 
follow for t > 0. We conclude that p = \{fl + aa), solves the truncated initial value problem for 
(2.1) globally in time. 
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We can remove the truncation using a variant of the particular limit procedure from [PLL] for 
Pj, when j — > oo. 



Theorem 2.2 [A2] Suppose that (/„, a ) € L°° n L^fO, l] 3 x M 3 ), and /or p 2 < f, j G N, i/iai 
< ?7j < /o < 2 — and a 2 + n 2 < min(/g, (2 - / ) 2 ) / or some rjj > 0. Taen iae system (2.2- 
3) with initial value (/o,cro) > a bounded integrable solution for t > wi/t < / < 2 and 
a 2 <min(/ 2 ,(2-/) 2 ) . 

For details see [A2]. Other problems of considerable physical interest are (cf [JM]) the relaxation 
times for spin-diffusion, the time asymptotic behaviour in general, and the influence of more involved 
transport terms in (2.2-3) such as the physicists' version of the problem, namely 



dt 



+ Vp £ p " Vr Jp Vr£p Vp 7p + [ Qp. ' Q r . 

df p dh p 



dh p da v 



— Qn 



de p dop 
dpi dri 



i=xyz 

depdap dfpdhp _ 
dri dpi dri dpi dpi dr { 



dn dp L 
2(h p x dp) = Q 



Here 



P 

2m 



+ 



J dp'{V(0)-±V(\p-p'\)}f p/ (r,t) 

h P = -^B- J dp'V(p-p')G p ,{r,t), 

V is the inter-particle potential, B an external magnetic field, and 7 is the gyromagnetic ratio. 

In spintronics for semiconductor hetero-structures, related linear spinor Boltzmann equations 
d 



dt 



p + v- \/ x p + E ■ \/ v p = Q(p) + Qso(p) + Qsf(p)- 



(2.7) 



are considered. Here E is an electric field and Q is the collision operator for collisions without 
spin-reversal, in the linear BGK approximation 



/ 



a{v,v')(M(v)p{v') - M(v')p{v)dv', 



M denoting a normalized Maxwellian. The spin-orbit coupling generates an effective field f2 making 
the spins precess. The corresponding spin-orbit interaction term Qso{p) is given by ■ <J,p]. 
Finally Qsf(p) is a spin-flip collision operator, in relaxation time approximation given by 



Qsf(p) 



tr pl 2 - 2p 



with T s f > the spin relaxation time. 

Mathematical properties such as existence, uniqueness, and asymptotic behaviour, have been stud- 
ied in particular by the French group around Ben Abdallah, his coworkers and students (see [BH] 
and references therein). 
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3 a) Low temperature bosons and a condensate. 



For a Bose fluid and for temperatures sufficiently low for the mass of the condensate noticeably to ex- 
ceed the mass of the (quasi-)particles, only interactions between thermally excited (quasi-)particles 
and a condensate are of practical importance. We shall now consider one such situation allowing 
transfer of atoms between the components, so that in particular no mass conservation should be 
expected for an individual component. The model (cf [KK], [K], [E] and for validation [SI]) is based 
on the Beliaev-Popov approximation of the Bogoliubov Green-function description which ignores 
off-diagonal correlations, and the Thomas-Fermi approximation which neglects a quantum pressure 
term. 

Our simplified two-component model consists of a kinetic equation for the distribution function 
of a gas of bosonic (quasi-)particles interacting with a Bose condensate, which in turn is 
described by a Gross-Pitaevskii equation (cf [PS]). 
In the local rest frame, the kinetic equation for the excitations is 

^ + i v P (E P + v s p) • Vx/ ~^Vx (E p + v s p) ■ sj p f = ^Q(f, n c ). (3.1) 

Here v s is the superfluid velocity, E p the Bogoliubov excitation energy, and e a small parameter. 
The collision term becomes 

Q(f,n c )(p) = n c J \A\ 2 5( Pl -p 2 - P3 )5(E 1 -E 2 -E 3 )[S(p- Pl ) (3.2) 

-5(p- P2 ) - 5(p-p 3 )]((l + h)hh ~ /i(l + / 2 )(1 + h))d Pl d P2 dp 3 . 

We notice that the domain of integration is only 2d and not 5d as for classical Boltzmann. The 
transition probability kernel |yl| 2 can be explicitly computed by the Bogoliubov approximation 
scheme. 

A := (u 3 - v 3 )(uiu 2 + viv 2 ) + (u 2 - v 2 )(uiu 3 + viv 3 ) - (u\ - v{)(u 2 v 3 + v 2 u 3 ). 
Here the Bose coherence factors u and v are 

2 £p + E p 2 2 i 

u p = ~^Ep~> u p- v p = 1. 

2 

with lp = ^ + gn c , n c the non-equilibrium density of the atoms in the condensate, m the atomic 
mass, and g a scattering length defined later. 

The collision operator Q(f,n c ) can be formally obtained (cf [ST], [EMV], [No]) from the Nordheim- 
Uehling-Uhlenbeck collision operator 

QnuuUXp) = I B5( P + p* - (p' + p',))5(E(jp) + Efa) - (E(jp') + E&))) 

JR 3 xl 3 xl 3 

(/7:u + /)(i + /*) - //*u + /0(i + fi))dp,d P f d P :. 

Namely, assume that a condensate appears below the Bose-Einstein condensation temperature T c . 
That splits the quantum gas distribution function into a condensate part n c 5 p= o and an L 1 -density 
part f(t,x,p), and gives 

QnuuU + n c S p= o) = Qnuu(I) + Q(f, n c ) + bn 2 c + cn 3 c + dn c S p=0 , 



7 



where a simple computation shows that b = c = 0. 

At the low temperatures we have in mind here, the number of excited (quasi-)particles are consid- 
ered to be small, when they are sufficiently excited for pair collisions to be important. Also the 
time-scale to reach equilibrium for such collisions, is considered to be short compared to the one for 
Q, and so the Qnuu collision term is neglected relative to the collision operator Q. 

The usual Gross-Pitaevskii (GP) equation for the wave function ip (the order parameter) 
associated with a Bose condensate is 

dib h 2 

ih ~it = ~2m~ Ax ^ + ( Uext + 5| ^ 2) ^' 

i.e. a Schrodinger equation complemented by a non-linear term accounting for two-body inter- 
actions. U ex t is an external potential. For simplicity, we here do not include the strongly space- 
inhomogeneous trapping potential, otherwise an essential ingredient in experiments on laser-trapped 
Bose gases. Modulo a numerical factor, g is the s-scattering length of the two-body interaction po- 
tential. In the present context the GP equation is further generalized by letting the condensate 
move in a self-consistent (Hartree-Fock) mean field 2gh = 2g J f(p)dp produced by the thermally 
excited atoms, together with a dissipative coupling term associated with the collisions. This gives 

i^t = -±A x i; + (U ext +g\i;\ 2 -i^ j C(f,n c )(p)dp)i>, (3.3) 

with C(f, n c ) = ^Q(f, n c ) and U ex t in our case equal to a constant plus a term of order e. It is often 
useful to split the equation for tp = ^/n^e 10 into phase and amplitude variables (polar representation 
or the Madelung transform) and remove some physically negligible term, leading to 



e2 -gf + Vx • {n c v s ) = - Q(f, n c )dp 



86 mv 2 



^7 = -A»c 



dt ^ c 2e 2 ' 
with fi c a local condensate chemical potential. 



3 b) A space-homogeneous, isotropic case. 

Consider the Cauchy problem for the two component model (3.1), (3.3) in a space-homogeneous 

isotropic situation and in the superfluid rest frame (condensate velocity v s = XJ X 9 = 0), i.e. the 
equations 

^ = Q(f,n c ), (3.4) 

^ = ~J Q(f,n c )dp, (3.5) 
with initial values 

f(p,0) = fi(\p\), n c (0) = n ci . (3.6) 

Here f(p,t) is the density of the quasi-particles, n c {t) the mass of the condensate, and the collision 
operator Q is given by (3.2). 
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There are two physically important regimes (cf [E]). One is the very low temperature situation 
with all \pi\ « po (T < 0.4T C in the set-up of [E]), i.e. where physically all quasi-particle momenta 
are much smaller than the characteristic momentum po = \J2mgn c for the crossover between the 
linear and quadratic parts of the Bogoliubov excitation energy of the quasi-particles; 



E <?> : = i/^+^-w + s^rj = «wd + (3.7) 

Here c := yj ^ the speed of Bogoliubov sound. Setting m = gives po = c. In applications with 

\p\ << Po> the right hand side of (3.7) is usually taken as the value of E(p). 
The Bose coherence factors can then be taken as 



gn c 1 E(p) I gn c 1 E(p) 



2E(p) 2 V 2#n c ' p V 2£(p) 2 V 2ffra, 



f 9 \/ 9™' ' Up ~~ A/ 9Rf.l 9 \/ 9™' ' U P V P ~ lj 



which gives 



1 ^E(p,)E(p')E(pj) f^T I EjpQ I E{p>) / E(pJ 

2§ ( 5 n c )§ V 2 l V E(p*)E(p>) + V E(p*)E(ti) \J E{p')E{^)>- 

And so recalling that = -E(p') + -E(K), we obtain 



22 (5«c) 2 
With this A, the collision operator becomes 

Q(f,n c )(p)= [ x E{pi)E ^ E{p3) S( Pl -p 2 -p 3 )5(E 1 -E 2 - E 3 )[5(p- Pl ) 

J g 6 K 

S(p - P2) - 5{p - ps)]((l + h)hh - /i(l + / 2 )(1 + h))d Pl dp 2 dp 3 , (3.8) 
where % denotes the truncation for < A, 1 < i < 3 for a given A > 0. 

The opposite limit where all momenta \pi\ » po has the dominant excitation of (Hartree-Fock) 
single particle type (in [E] for moderately low temperatures around 0.7T C ). Expanding the square 

2 

root definition of E in (3.7), we may approximate E p by + gn c leading to a collision operator of 
the type 

Q(f,n c ) = kn c / xKPi -P2~P3)S(E 1 - E 2 - E 3 )[8{p-pi) 

Jr 3 xR 3 xR 3 

Sip ~ P2) - S(p - ps)]((l + h)hh - /i(l + / 2 )(1 + h))dpidp 2 dpz (3.9) 

for the 'partial local equilibrium regime'. Here x is the characteristic function of the set of 
ip,Pi,p 2 ,p 3 ) with \pi\, \p 2 \, \p 3 \ > a for a given positive constant a. 

In the general case, the collision operator is 

Q(f,n c )ip) = n c j \A\ 2 5( Pl - p 2 - p 3 )5(E 1 - E 2 - E 3 )[8(p - Pl ) 
Sip - p 2 ) - Sip - p3)]((l + h)f 2 f 3 -hil + f 2 )il + h))d Pl dp 2 dp 3 , (3.10) 
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with the excitation energy E defined by 
The kernel |vl| 2 is bounded by a multiple of 




in the physically interesting cases where asymptotically all \pi\ « po, all \pi\ >> po, or one 
\pi\ « po and the others >> pq. These three cases are relevant for very low respectively moderately 
low temperatures compared to T c , and (the third case) for the collision of low energy phonons with 
high energy excitations (atoms). The asymptotic situation of two \p-i\ « po and one pi >> po 
(with unbounded A) is excluded by the energy condition. 

Using \A\ 2 as the kernel in the collision operator, we have proved the following result. 

Theorem 3.1 [AN1] Let n c i > and fi(p) = fi{\p\) G L 1 be given with f\ nonnegative and 
fi(p)\p\ 2+,y £ L 1 for some 7 > 0. For the collision operator (3.10) with the transition probabil- 
ity kernel \A\ 2 , there exists a nonnegative solution (/, n c ) G C 1 ([0, 00); L l + ) x C 1 ([0, 00)) to the 
initial value problem (3.1-3) in the space- homogeneous, isotropic case. The condensate density n c 
is locally bounded away from zero for t > 0. The excitation density f has energy locally bounded in 
time. Total mass Mq = n c i + J fi(p)dp is conserved, and the moment f \p\ 2+1 f'dp is locally bounded 
in time. In the moderately low temperature case a total energy quantity is conserved. 

The proof is via approximations controlled by a priori estimates and fixed point techniques. An 
important still open problem is the question of time-asymptotics. The addition of a NUU collision 
term would (using methods we are aware of) considerably weaken the type of solutions that can be 
obtained. 



A. Heintz has developed numerical methods for the above evolution problem in connection with 
an exam project [BDJ]. The recent paper [S2] considers the spatially homogeneous and isotropic 
kinetic regime of weakly interacting bosons with s-wave scattering. It has a focus on post-nucleation 
self-similar solutions. Another recent paper, [EPV], studies linearized space homogeneous kinetic 
problems in settings related to the ones discussed here, and with a focus on large time behaviour. 



3 c) A space-dependent, close to equilibrium case. 

In equilibrium, the right hand side of the kinetic equation (3.1), i.e. the collision term, vanishes. 
Multiplying with jq^^j and integrating in p, we obtain 



(1 + /*) (! + /')(! + /J) 
As usual - and ultimately because of a Cauchy equation - equation (3.11) implies that / = 

2 

with M a Maxwellian. We assume that the Maxwellian has zero bulk velocity, M = exp(— + 5). 
With pq = y/2mgn c , we restrict to the 'intermediate temperature region' \p\ >> po in the super- 
fluid rest frame. Denoting the mass density of the superfluid by no, and using the approximation 
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2~ + 9 n o for the energy E of (3.7), equation (3.11) holds for 5 = — -^-griQ. Hence for m = ^, the 
particular Maxwellian Ai$ = exp(— f3(p 2 + gno)) gives an equilibrium for the collision operator Q. 
Take for simplicity g = (3 = 1, which fixes the equilibrium limit as (Pq, uq) with Pq = ( cxp (pi+ no )_i) > 
the Planckian corresponding to A^o- Require that the ingoing boundary values for the kinetic part 
are equal to this limit. Take the initial values (fo,no) as e 3 -small deviations from this equilibrium 
which retain the total mass Mq. 



This space-dependent, close to equilibrium evolutionary problem poses interesting questions about 
existence (work in progress with A. Nouri) and the long time behaviour, including whether the 
solution converges to equilibrium when time tends to infinity. A collision term of NUU type can 
easily be included. In contrast to the NUU case of [R], here the splitting of the linearized collision 
operator into collision frequency plus a compact operator, requires a (physically acceptable) cut-off 
for very large velocities. Also the Milne problem is different from the classical case, e.g. in not 
having a constant mass flow, and again requires a truncation for large p's. 



4 a) A kinetic anyon equation. 

The quantum Boltzmann Haldane equation [BH] is a kinetic equation of Boltzmann type for confined 
quasi-particles with Haldane statistics [H]. This exclusion statistics interpolates between the Fermi 
and Bose quantum behaviours. In quantum statistical mechanics, the number of quantum states of 
TV identical particles occupying G states is given by 

(G + N-lV. , G\ 
and 



N\{G-l)\ N\(G-N)\ 

in the boson resp. fermion cases. The interpolated number of quantum states for the fractional 
exclusion of Haldane and Wu is 

(G+(JV -Dd-o)), o<a<i 



N\(G - aN - (1 - a))! 



When applied to the fractional quantum Hall effect, the Haldane statistics coincides with the two- 
dimensional anyon definition (in terms of the braiding of particle trajectories. 
Haldane statistics may also be realized for neutral fermionic atoms at ultra-low temperature in three 
dimensions. 



Elastic pair collisions in a Boltzmann type collision operator, are expected to preserve mass, lin- 
ear momentum, and energy. That holds for the collision operator Q of Haldane statistics, first 
introduced in [BBM], 

Q(f)= I B(v-v*,u;)xlf'flF(f)F(f*)-ff*F(f')F(fl)}dv*du. (4.2) 

JM d xS d - 1 

Here duo corresponds to the Lebesgue probability measure on the sphere. 

The collision kernel B in the variables (z,u) £ M d x S 6 '" 1 is positive, locally integrable, and only 
depends on \z\ and |(z,oj)|. The filling factor F is given by F(f) = (1 — af) a (l + (1 — a)/) 1- ", 
< a < 1. The factor (1 — af) keeps the value of / between and ^. F is convex with maximum 
value one at / = for a > ^, and maximum value (^ — l) 1 " 2 " > 1 at / = ^fj^ for a < \. 
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With this filling factor, the collision operator vanishes identically for the Haldane equilibrium dis- 
tribution functions as obtained in [W] under (4.1), but for no other functions. For the limiting 
cases, representing boson statistics (a = 0) and fermion statistics (a = 1) the quartic terms in the 
collision integral are canceled, and ideas from Lions' compactness result for the classical gain term 
are applicable. This allows approaches related to the weak L l analysis for the quadratic Boltzmann 
equation. 

But for < a < 1 there is no cancelations in the collision term. Moreover, the Lipschitz continuity 
of the collision term in the Fermi-Dirac case, is now replaced by a weaker Holder continuity. For 
the space homogeneous case, I have developed an approach based on strong L -'--compactness. 

Consider the initial value problem for the Boltzmann equation with Haldane statistics in the space- 
homogeneous case with velocities in M d , d > 2, 

| = Q(/). (4.3) 

Because of the filling factor F, the range for the initial value /o should belong to [0, which is 
then formally preserved by the equation. The general BH equation (0 < a < 1) retains important 
properties from the Fermi-Dirac case, but it has so far not been validated from basic quantum theory. 
In particular the choice of kernel for the collision operator has to be studied in other ways. The choice 
by the physicists has been to use classical kernels such as those for hard forces, B(z, oj) = |z| /3 6(^j^) 
with < (3 < 1 and Grad angular cut-off, which I will now discuss. That also agrees with the better 
understood limiting cases a = 0, 1. We assume that < b < c\z\@\ sin 6 cos with an initial 

moment (1 + \v \ s )fo in L°° for s > d — 1 + j3. 



Theorem 4.1 [A3] Consider the space-homogeneous equation (4-3) for hard force kernels with 

< B(z,u) < C\z\P\ sm6 cos e^ 1 , (4.4) 

where < /? < 1, d > 2, and < /? < 1, d = 2, and let the initial value < /o € L 1 have 
finite energy. If < /o < - and esssup(l + |w| s )/o < 00 for s = d — 1 + ft, then the initial value 
problem for (4-3) has a solution in the space of functions continuous from t > into L 1 n L°° , 
which conserves mass and energy, and for to > given, has esssup„ t<to \v \ s f(t,v) bounded,where 
s ' = min(S) W+±2±E). 

Idea of proof. This initial value problem is first considered for a family of approximations with 
bounded support for the kernel B, when < /o < esssup/o < ^. Starting from approximations with 
Lipschitz continuous filling factor, the corresponding solutions are shown to stay away uniformly 
from ^, the upper bound for the range. Uniform Lipschitz continuity follows for the approximating 
operators and leads to well posedness for the limiting problem. Uniform L°° moment bounds for the 
approximate solutions hold, using an approach from the classical Boltzmann case [Al]. Based on 
those preliminary results the global existence result for hard forces of Theorem 4.1, can be proved 
by strong compactness arguments. Mass and first moments are conserved and energy is bounded 
by its initial value. That bound on energy in turn implies energy conservation using the arguments 
for energy conservation from [MW] or [Lu2]. ■ 
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Remarks. 

i) The proof of the theorem implies the following stability. Given a sequence of initial values (/on)neN 
with 

supesssup v fon(v) < — , 
n a 

and converging in L 1 to /o, then there is a subsequence of the solutions converging in L 1 to a 
solution with initial value Jq. 

ii) The question of long time behaviour is open. 

hi) The space-dependent case is of interest also in the perspective of the general existence theory 
for Boltzmann type equations. 



5 Differences between quantum and classical Boltzmann theory. 

The examples discussed above , were selected to illustrate how the situation in quantum kinetic 
theory in various ways markedly differ from classical kinetic theory. The selection could be ex- 
panded to other types of low temperature kinetic evolutions of collision type, such as collisions 
involving five quasi-particles, two colliding ones giving rise to three or conversely, cf [K]. Other 
examples are Grassmann algebra valued gas densities, i.e. completely anti-commuting densities, 
e.g. for modelling the Pauli exclusion effect, cf [PR], and pair collision terms taking account of the 
finite duration of a collision, cf [ES]. 

i) Qualitatively different collision operators. 

The type of collision operator varies qualitatively much more in the quantum regime than in the clas- 
sical one, like (1.1), (1-2), (2.1), (2.7), (3.2), (4.2), and the examples mentioned from [K], [PR], and 
[ES]. The quantum influence sometimes appears as entirely new phenomena such as the collision- 
driven spin waves in the fermionic spin matrix Boltzmann equation. The thermal de Broglie wave 
length may become much larger than the typical inter-particle distance, an aspect totally absent 
from the classical theory. The quantum influence leads to different filling factors which, as for the 
NUU equation, can stabilize as well as destabilize the collision effects. 

ii) Fewer collisions. 

There are fewer collisions close to absolute zero, where permitted energy levels usually are discrete. 
In the example of boson excitations plus a condensate, an excitation may interact only with the 
modes of the zero-point motion that do not give away energy to it. Also, in that example the 
domain of integration has lower dimension than the classical case. 

iii) Typical energy range. 

A quantum situation is usually not scale invariant as in the classical case, but may have a typical 
energy range (Section 3b). A particular type of quasi-particles/excitations exists in a particular 
energy interval, and extending the energy outside this interval, may introduce collision types not 
observed in the experiments, and physically irrelevant for the modelling at hand. In this way 
bounded velocity domains may be both physically correct and mathematically adequate in quan- 
tum situations, like the two component, space-dependent boson example. Similarly, models and 
solutions with respect to a bounded interval of time, e.g. depending on relaxation effects, may be 
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relevant in particular quantum situations, such as the exciton-polariton case of [DHY]. 

iv) Parameter range and kernel family. 

The collision terms and energy expressions discussed, were defined from a physical context. Chang- 
ing a parameter mildly may completely change the mathematical aspects of the problem, as in the 
anyon example for a around and 1. Another example where very small changes in a parameter 
completely changes the physical situation, is the the extremely narrow temperature-pressure domain 
for the phenomenologically very rich A-phase in 3 He. The question of parameter range and kernel 
family, is in general more delicate for quantum than for classical kinetic theory. 

v) Questions from classical kinetic theory. 

For each type of quantum kinetic collision term, the multitude of questions from classical kinetic 
theory can be posed and studied. New boundary-condition dependent phenomena may appear, 
an example being a low-temperature gas of excitations in a condensate between rotating cylinders 
having a vacuum friction type of radiation as in the Zeldovich-Starobinsky effect, cf [V]. 

vi) New insights. 

The quantum problems may - in comparison with their classical counterparts - require new ap- 
proaches or additional ideas for their solution, like the Milne problem for the space-dependent boson 
system. 
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